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ABSTRACT 

Aims. It is commonly assumed that the two indirect exoplanet detection methods, the radial velocity method and as- 
trometric method, require observational periods exceeding the orbital period to produce positive results. Here we test 
this assumption in detail. We also investigate the smallest ratio of observational timeline and orbital period required 
for positive detections. 

Methods. We obtain full information on the orbital parameters by combining radial-velocity and astrometric measure- 
ments by means of Bayesian inference, and sample the parameter probability densities of orbital and other model 
parameters with a Markov Chain Monte Carlo (MCMC) method in simulated observational scenarios to test the de- 
tectability of planets with orbital periods longer than the observational timelines. 

Results. We show that, when fitting model parameters simultaneously to measurements from both sources, it is possible 
to extract much more information from the measurements than when using either source alone. Currently available 
high- precision measurements of radial velocity (with 1ms -1 precision) and astrometric measurements achievable with 
the SIM space telescope (with a precision of lfias) can be used together to detect a Jupiter analog 30pc away with an ob- 
servational timeline of only three years, approximately one fourth of the orbital period. Such measurements are sufficient 
for determining all its orbital parameters, including inclination and the true mass. Also, with accurate radial velocity 
measurements covering a timeline of 20 years, the true mass could be determined by astrometric observations within a 
single year. These case studies demonstrate the potential power of the Bayesian inference of multiple data sources in 
exoplanet observations. As an example, we show that using the currently available radial velocity measurements, the 
inclination of HD 154345b could be determined with SIM in a year. 

Key words, planetary systems - astrometry - Methods: statistical - Techniques: radial velocities - Stars: individual: 
HD 154345 



. ^ 1. Introduction to detect the inclination, and consequently the true mass 

^^ . of the companion. But despite several trials (van de Kamp 

"£ . Since the discovery of 51 Peg b (Mayor & Queloz 1995), 1969; Neuhaeuser et al. 2008), the precision of these mea- 

C3 ; high-precision radial velocity (RV) measurements have been sureme nts has not been high enough to verify the planetary 

used successfully to detect planetary companions of nearby nature of RV companions 

stars (Butler etal. 2006). The finest instruments to date are with thc aid of telescopes and instruments capable 

capable of achieving an RV precision of 1ms (Santos et al. of h i gn _ prec i s ion astrometry (SIM, GAIA, PRIMA, etc.); 

2004; Moorwood and Masanon, 2004). However, the exact howeverj this situ ation is about to change. With an esti- 

nature of the companions remains unknown when only RV mated astr0 metric precision of future telescopes of l.O^as 

measurements are available, as these only yield the prod- (Unwin et al 2Q08) or 8 _ 10 ^ as (Casertano et al. 2008; Derie 

uct of mass and sine of the orbital inclination, giving the et aL 2002); it wiU be possible to determine the masses of 

lower limit for the mass. Recently, it has been claimed that the already detected RV companions. It is commonly as- 

an RV precision of lems could be technically possible m sumed that these det ections require observational periods 

the future (Li et al. 2008). However, the RV variations of longer than the orbital period of the target system to be 

stars of approximately 5ms to 50ms for K to F stars, able to detect the periodic signaL Here we test this assump- 

caused by star spots and irregular convective zones, will ^ on j n de t ad 

likely prevent the detection of signals of the cm scale (Saar " In tMg artide we simulate astrom etric and RV measure- 
& Donahue 1997; Saar et al. 1998). ments to gtudy the possibility of detecting planetary com- 
Astrometric (A) measurements of the position of the pan i on s of nearby stars with various observational time- 
target star in thc sky as a function of time can be used lines. T h e goal is to find the minimum timeline required 

for detecting a planetary companion using high-precision 

* e-mail: mikko.tuomi@utu.fi RV and astrometric measurements. These two sources of 

** e-mail: samuli.kotiranta@utu.fi data are combined by means of Bayesian inference, and the 
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probability densities of orbital and inertial reference frame 
parameters are sampled using Markov Chain Monte Carlo 
(MCMC) (Metropolis et al. 1953; Hastings 1970) to find 
the full global solution to this multidata inverse problem 
(see Ford 2005; Ford 2006; Maness et al. 2007). Also, the 
probability densities are sampled to calculate the realistic 
error bars for parameters and to determine whether a pos- 
itive detection has been made or not. We also analyse the 
RV measurements of HD 154345 (Wright et al. 2008) alone 
and together with simulated astrometric measurements to 
estimate the minimum observational timeline of astrometry 
for the detection of the true mass of HD 154345 b. 



2. Modelling the data 

The motion of a planet around a star was treated as a sim- 
ple two-body system, with masses m* and m p for the star 
and the planetary companion, respectively. In Cartesian co- 
ordinates, when the gravitational forces between the pos- 
sible other planets in the system are assumed negligible, 
the column position vector of the star with respect to the 
barycentre of the system can be expressed as a function of 
time (t) as (e.g. Green 1985) 



R(t) = R(Q) + R(0 )t + P[cos E{t) - e] 
+QVl~e 2 sin E(t), 



(1) 



where E is the eccentric anomaly, satisfying the Kepler 
equation E — esvtxE — 2irn{t — to), and n = P^ 1 is the 
orbital frequency, P the orbital period and to the time of 
periastron. The velocity and position R(0) and i?(0) w.r.t. 
the observer are some constant vectors defining the iner- 
tial reference frame. Vectors P and Q are constant vectors, 
defined as 



P = a*(i cos u> + k sin u>) 
Q = a*(— I sin a; + k cos w) 
I = (sin.0, cosO, 0) 
k = (I cos il,—I sin f2, sin i) 



(2) 



The parameters a*, e, Cl, I = cos«, and u> are the orbital 
parameters of the system: semimajor axis of the star, eccen- 
tricity, longitude of ascending node, cosine of the inclination 
and the longitude of pericentre, respectively. The semima- 
jor axis of the star can be expressed as a function of the 
masses of the gravitationally interacting bodies, 



a* = mpG 1 / 3 ^™* + m p \) 2/3 . 



(3) 



The models for astrometric and RV data are now sim- 
ply &(t) = [Q x (t),Q y (t)} = D- l [R x {t),R y {t)] for as- 
trometry and z{t) — dtR z (t) for RV. Here D is the dis- 
tance of the system from the observer. There are now 12 
independent parameters describing the system. The pa- 
rameter vector u of parameter space U can be written 
as u = (X x ,Xy,/j, x ,fj, y ,'Y,a it ,LJ,I,0,,e,to,n) G U, where 
A, = D- l R x (0), X y = D- l Ry(0), fx x = D^R^O), /i y = 

D~ 1 R y (0) and 7 = R z (0). To fully describe the system, 
the masses m p and m* should be treated as independent 
unknown parameters, and the distance D should also be 
included in vector u. Here we assume that m+ and D are 
known with sufficient accuracy by some astrophysical tech- 
niques. 

It is essential to include the parameters defining the 
inertial reference frame to be able to fully investigate the 



probability densities of the parameters. It is also assumed 
that the model parameters do not change as a function of 
time during the observational timeline. 

It is intuitively clear that the availability of both RV and 
A data sets, instead of only one of them, should increase 
the amount and the quality of information on the observed 
system. This should be true even if one of the data sets is 
significantly more inaccurate than the other. In what fol- 
lows, we state this principle more rigorously and present 
a practical method for combining the two data sources. In 
fact, since RV and A data are radial and tangential projec- 
tion samples of orbital motion, they are separate sources 
only in the instrumental sense, so we could just as easily 
talk about full velocity data. 



3. Simulations and fitting procedure 

We simulated data sets to study the complementarity of as- 
trometric and RV time series. Assuming that the errors ej 
for each data point i = 1, ..., N are independent and iden- 
tically distributed and that their probability density func- 
tions (PDF's) are Gaussian (a ~ N(0, of)), the observed 
RV data are modelled as 



Z(u,ti) + £j. 



(4) 



Equally, the two-dimensional astrometric model can be 
written as 



0, =&(u,U) 



(5) 



To avoid inversion crimes, i.e. to avoid using exactly the 
same model to generate the measurements and to find the 
inverse solution (see e.g. Kaipio and Somersalo 2005), an 
additional planet of low mass was included in the model 
when generating the data. This choice was made because 
small systematic errors make the simulated measurements 
more realistic. Otherwise the simulated measurements and 
the corresponding solutions would only help for studying 
the model in Eq. ([1]) , not necessarily situations encountered 
in reality. 

Several data sets were generated, each with a different 
value for the data parameter T, representing the length of 
the observational timeline. It was further assumed that the 
observation time ti was evenly distributed in T. 

When simulating the measurements, the timelines were 
set to T A = T RV = 20.0, 10.0, 5.0, 3.0, 1.5, 1.2, 1.0 and 0.8 
years. The values of all the other data and orbital param- 
eters and masses were fixed in all these scenarios. These 
values were set to (<jkv , &a, Nrv , Na) — (l.Om/s, 1.0/ias, 
100, 100) and (a,e,u>,i,Cl,to,m p ,m i ,,D) = (5.0AU, 0.1, 

1.0, 1.0, 1.0, lOOO.Od, Mj, M Q , 30pc), where Mj is the 
mass of Jupiter. This simulated system is called SI. 

3.1. What is a positive detection? 

In the simplest possible case, when e = Oand i — 7r/2, the 
detection threshold of full velocity data can be calculated 
analytically. Let N be the number of data points, T the 
length of the timeline of observations, and a the standard 
deviation of observations. Following the approach in Eisner 
& Kulkarni (2001a, 2001b), the detected signal of the ve- 
locity variation amplitude is a false one produced by the 
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uncertainties in the data with a probability < 1% if 




(G) 

(7) 

(8) 

(9) 

(10) 



and a a is in radians. This approach excludes the uncer- 
tainties in the orbital period and can therefore only yield 
the lower limit for the detection threshold. Hence, if Eq. 
([6]) does not hold, it will be impossible to detect the signal. 
However, if it holds, the detectability of such a companion 
needs to be examined more closely by numerical simula- 
tions and by analysing the simulated data using methods 
such as MCMC and Bayesian model selection criterion. 

To fully investigate the ability to detect planetary com- 
panions, we must define when a positive detection has been 
made. This question can be approached through Bayesian 
probabilities. Let i?i be the model in Eq. (Q]) with one plan- 
etary companion (corresponding 12 parameters in the RV 
and astrometry models), and Rq a model without a plan- 
etary companion (5 parameters). In general, let R k be a 
model with fc planets. 

Using the Bayes theorem, it can be seen that the con- 
ditional probability of model Rk representing the data (m) 
best, out of the p + 1 alternatives to be tested, can be writ- 
ten as 



P(R k \m) = P(R k 



J2B k Am)P(Rj) 



i=o 



(11) 



where the Bayes factor B k j 
Raftery, 1995) 



Bfcj(m) 



P(m\R k ) 
P(m|H,-)' 



is defined as (e.g. Kass and 



(12) 



and P(R k ) is the prior probability of the fcth model, here 
set equal for all fc, because it is assumed that there is no 
prior information available. Here the likelihood P(m\R k ), 
with parameters u k g U k for the fcth model, is 



P(rn\R k ) 



UktUk 



p(m\u k , R k )p(u k \R k )dUj , 



(13) 



where p(m\u kl R k ) is the parameter likelihood function and 
p(u k \R k ) the prior density. 

Since the model probability, defined in this way, auto- 
matically takes the Occamian principle of parsimony into 



account, the model with the smallest number of parame- 
ters out of those having almost equal probabilities will be 
selected. Hence, it can be said that a detection has been 
made if (Jeffreys 1961) 



P( J Ri|m)>P( J R |m). 



(14) 



This criterion is used throughout this article when deciding 
whether a statistically significant detection has been made 
or not. 



3.2. Fitting method 

The fitting was performed by requiring that the values of 
all the least-squares cost-functions S x (astrometric x), S y 
(astrometric y), Srv (RV), and their sum be minimized si- 
multaneously. This method, called multidata inversion, has 
been used successfully with astrometric and RV measure- 
ments when detecting stellar binaries (e.g. Torres 2007). See 
the discussion in Kaasalainen and Lamberg (2006), where 
the multidata inversion was applied to asteroid observa- 
tions. 

The models for astrometric position and RV of the 
two-body system of interest are non-linear, so an itera- 
tive method of fitting the model parameters is needed. 
The MCMC with Metropolis-Hastings (M-H) algorithm 
was chosen because it is a global method (Metropolis et 
al. 1953; Hastings 1970), it offers a direct estimate of the 
posterior probability density, and because it can be used to 
verify the existence and uniqueness of the solution. Since 
the probability densities given the measurements are avail- 
able, MCMC can be used to calculate realistical error esti- 
mates for the model parameters. These estimates are typ- 
ically much larger than those calculated using traditional 
methods (e.g. Ford 2006), implying that MCMC should be 
preferred when assessing the parameter errors. Assuming 
Gaussian errors with zero mean, the likelihood function of 
the parameters with respect to RV measurements can be 
written as 



p(z\u) ex exp(— 0.5Srv)- 



(15) 



When applying MCMC, a parameter value (tto) is selected 
for the first member of the chain. The next value Uk+i 
is found by randomly selecting a proposal in the vicinity 
of Uk ■ This is then accepted by comparing the likelihoods 
of the two parameter values. Proposed parameter values 
Wfc+iwith a greater likelihood than that of u k are always se- 
lected as the next chain member, but values with a smaller 
likelihood can also be selected according to the criterion of 
Hastings (1970). Samples of at least 10 5 points were gen- 
erated when sampling the parameter space. For practical 
details on MCMC with astronomical data, see e.g. Gregory 
(2005). 

The parameter space U in this Keplerian two-body 
model has a comparatively small dimension (dimU = 12), 
but in some cases it already makes the sampling compu- 
tationally expensive. Especially when covariances between 
the parameters are large and of non-linear nature, the space 
of reasonable probability Ur C U to be sampled can be very 
narrow and, as a result, the next proposed value of param- 
eter vector u in the Markov chain is likely to be outside 
this subspace and thus rejected, considerably increasing the 
time needed to generate a statistically representative chain. 
For this reason, when using a multivariate Gaussian density 
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as a proposal, the acceptance rates were low, approximately 
0.1 in the MCMC samplings. 

4. Taking advantage of Bayesian inference 

With more than one source of measurements available, it is 
possible to get more information from the system of interest 
than when relying on any single observation method alone. 
This is a consequence of Bayesian inference. 

Denoting the astrometric measurements by O and the 
RV measurements by z , the conditional probability of hav- 
ing parameter vector u G U, is a product of the impacts 
the two sets of measurements have on this hypothesis: 



p(u\z ,B ) = —— — L-ip( tt ). 

p(z ) p(@ ) 



(16) 



Thus, there is always more information available on the sys- 
tem - either in a narrower parameter density or in the pos- 
sibility of including more parameters in the model - when 
using multiple data sources. This is due to complementary 
rather than just additional information: the separate prob- 
ability densities for complementary sources are quite differ- 
ent from each other, so their product (joint probability) is 
much more tightly bound than either factor alone. 

Equation (|16p is in fact just another way of stating that 
we simply minimized the sum S = S x + S y + Srv, as can be 
seen by applying the Eq. (fT5j) . while checking that each Si 
was still close to their minima. However, when calculating 
model probabilities and parameter densities, the formula- 
tion in Eq. (fTl)]) has to be used. 

4.1. Correlations and complementarity 

By having measurements made using two different obser- 
vational techniques has an effect on the parameter PDF's. 
This happens because the two measurements are modelled 
using a different model with differing parameters for the 
inertial reference frame. Therefore, it is expected that the 
two measurements contain complementary information on 
different aspects of the system. 

Generally, correlations between model parameters occur 
if, for some small displacement of the parameter vector, Su, 
the model R used to describe the system satisfies 



R(u) - R(u + Su) w 0. 



(17) 



Now the parameter PDF's are broadened or correlated until 
Eq. (fT?) no longer holds. 

In a stellar system with a single planetary companion, 
the most obvious possible coupling, a positive correlation 
between to and u), is a natural byproduct of the two-body 
Keplerian model. When assuming e ~ and using Eq. |J3J), 
Eq. (fT]) becomes 



R(t) ~ m[l cos(cj + nt — nto) 

+k sm(uj +nt- nto)] + R(0)t + R(0). 



(18) 



Setting R(Lu,to)-~R(u>+5u,to + 5to) = implies 6u> = nSto, 
resulting in a positive linear correlation between parameters 
u> and to- 

If it is also assumed that the observational timeline is 
much shorter than the orbital period, T <C P, more cor- 
relations take place in this long-period system. As U £ 
[-T/2,T/2] for all i = 1,...,JV, when T/P -> 0, each 





a, [1e9 it,] 



trad] 



Fig. 1. Equiprobability contours containing 99%, 95%, 
90%, and 50% of parameter PDF's showing the densities 
of and correlations between parameters a* and n and pa- 
rameters / and f2. The simulated system has a long-period 
planet with T A /P = T RV /P « 1/4 (T RV = T A = 3.0 
years) . 



U/P — > as well. Thus, this assumption justifies cos(nt) 
1 and sin(ni) ~ nt and Eqs. ^ and ([5]) become 



&(t) « D- l M(P^ + Qfnt) 
z ~ Q$, z - P$, z nt + 7 



At 



M 



(19) 



where P$ and Q^ are just the functions presented in Eq. 
(|2j) with the angle to + nto replaced with 0, and matrix 
M = Af 3x2 = diag(l,l). 

Clearly, Eq. ([17]) holds if SP^ = —8/j, and nSQcf, + 
QcfiSn = —SX. This means that it is possible to change 
the values of the components of vector P$ by any amount 
and a corresponding negative change in the components of 
vector fi cancels this change exactly. As a result, the com- 
ponents of these vectors can correlate negatively. Also, the 
components of Q^ can correlate similarly with the com- 
ponents of A. For RV, Q^^z can correlate with 7, but the 
product Pcjt^n has no corresponding parameter to correlate 
with. These are just the correlations described by Eisner & 
Kulkarni (2001a, 2001b, 2002). 

From Eq. (fT9|) it is also clear that the orbital frequency 



can correlate with P$. z and Q,p making the detection of 
planetary signal harder and broadening the densities of the 
corresponding parameters. 

Despite the existence of and due to the partially comple- 
mentary nature of these correlations, it is possible to detect 
the periodic signal of a long-period planetary companion 
in the Bayesian model selection sense. The equiprobabil- 
ity contours of parameter combinations (n, a*) and [I, CI) 
with the simulated system SI are shown in Fig. [T] for 
Ta = Try = 3.0 years, which is approximately one fourth 
of the orbital period. The Bayesian model probabilities were 
found to satisfy the condition of Eq. [Mj and the contours in 
Fig.Q]demonstrate that it is indeed possible to detect extra- 
solar planetary companions even if the observational time- 
line is shorter than the orbital period. Also, since clearly 
I < 1, the planetary nature of these companions can be 
verified in this scenario. 

When using the two data sources, the planetary signal 
could not be detected for Ta — T R v < 3.0 years. For astro- 
metric or RV measurements alone, this signal was found to 
be undetectable for Ta = Try < 10-0 years, which clearly 
demonstrates the advantages of the Bayesian inference of 
multiple datasets. 
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Fig. 2. PDF's of parameters / and O in a snapshot scenario 
with T A /P w 1/11 (T A = 1.0 years). The mode, mean (/A 
standard deviation (ct), skewness (^ 3 ), and kurtosis (fi 4 ) 
of the densities are shown. The solid curve is a Gaussian 
function N([i,<t 2 ). 
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Fig. 3. Simulated RV and astrometric measurements and 
the maximum a posteriori orbit. 



4.2. Astrometric snapshots and detection thresholds 

Astrometric observations with the property Ta < P < Try 
are called astrometric snapshots. This definition is made 
because the astrometric observations are now made in a 
fraction of the time interval of the RV observations inca- 
pable of separating m and sini. We modified the simulated 
system SI by fixing Try — 20.0 years and denoting this by 
S2. 

Using Bayesian inference between RV and astrometric 
measurements made it possible to fit all the parameters in 
the model, including / and fi, to the measurements, even 
when the signal of the planetary companion could not be 
detected using astrometric observations alone. For the sim- 
ulated system S2, the condition in Eg . [T4l was satisfied for 
values of Ta as low as 1.0 years, which is less than one tenth 
of the orbital period. The corresponding densities of I and 
fl in this scenario are shown in Fig. [2j It can be seen that 
the maximum a posteriori estimates of these densities are 
very close to the values selected for the simulated system. 
This implies that the true mass of the stellar companion 
is obtainable. Also, the density of / shows that certainly 
/ < 1, implying that a companion of planetary mass has 
indeed been detected as claimed. The simulated RV and 
astrometric measurements and the maximum a posteriori 
orbit are shown in Fig. [3] for this snapshot scenario. 

The reason the parameters / and tt can be fitted is that 
the parameters in the RV model are now well-constrained 
by the RV measurements. The only possibility for Eq. (fT7|) 
to be true is that a+n sin i, the amplitude of RV variations, 
remains unaltered even if a* and i do not. This implies that 
a* ex (1 — I 2 ) -1 / 2 , which is the correlation observed in the 
parameter densities of a+ and / (Fig. 2J). Because of this 
correlation, the reference frame parameters of astrometry 
can correlate freely with a*. This is also demonstrated in 
Fig. SI 



Fig. 4. Equiprobability contours containing 99%, 95%, 
90%, and 50% of parameter probability densities. The cor- 
relations between a* and I and of a* and \ x in a snapshot 
scenario with T A /P « 1/1 1 (T A = 1.0 years). 




1000 



2000 3000 
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Fig. 5. RV measurements of HD 154345 and the maximum 
a posteriori orbit of the planetary companion. 



5. Case study: HD 154345b 

Recently, Wright et al. (2008, hereafter W08) reported a de- 
tection of a Jupiter analog orbiting a G8 dwarf HD 154345. 
They claim that astrometric measurements over its 9-year 
period would determine the orientation of the orbital plane 
and as a consequence the true mass. The data published in 
W08 was re-examined and the orbital solution found using 
MCMC. These data have 55 measurements over a period of 
10.4 years. The largest gap between two subsequent obser- 
vations within these measurements is 352 days. The orbital 
parameters were calculated assuming the same jitter level 
as in W08, 2.5ms" 1 , and are listed in Table [Q This Table 
shows the MAP estimates of the parameters and their 99% 
confidence sets (CS). Missing confidence sets indicate that 
the posterior density of the corresponding parameter has 
significant values everywhere in its parameter space. The 
results of Wright et al. (2008) with 99% confidence limits 
are shown for comparison. The corresponding fit is shown in 
Fig. [5l The large uncertainties of parameters u> and to (their 
99% Bayesian confidence sets are equal to their parameter 
spaces) stem from these parameters being meaningless for 
circular orbits. 

Regardless of the fact that in W08 the RV movement 
of the host star had been subtracted from the data, the 
constant movement parameter 7 was fitted as well to be 
able to take its uncertainty and its effect on the orbital pa- 
rameters into account. The orbital parameters were found 
consistent with but their confidence intervals smaller than 
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Table 1. The solution and error estimates of the HD 154345 
system. 



Parameter 


MAP 


99% CS 


Wright et al. 


P [years] 


9.06 


[8.56, 9.72J 


9.15 ± 0.67 


e 


0.02 


[0.00, 0.17] 


0.044 ± 0.118 


u>[°] 


90 


- 


68 


t [JD] 


245000 


- 


2452830 ± 850 


m p smi [M Jup ] 


0.95 


[0.80, 1.10] 


0.947 ± 0.232 


a [AU] 


4.16 


[4.01, 4.37] 


4.19 ± 0.67 


7 [ms _1 ] 


0.01 


[-1.42, 1.42] 






those reported in W08. This difference in the confidence 
intervals takes place likely because the Bootstrap method 
was used to assess the parameter errors in W08 instead of 
direct sampling of the posterior density. 

Astrometric measurements with a a — 1/ias, Na = 100, 
and T A =10.0, 5.0, 3.0, 2.0, 1.5, 1.2, 1.0, and 0.8 years were 
generated to study the detectability of parameters / and Q. 
These measurements were generated assuming that there is 
a planetary companion with orbital parameters in Table [TJ 
and (/, f2) = (0.0,1.0 rad). This simulated data was used 
together with the real RV measurements published in W08 
to find the limiting Ta for which the true mass of the planet 
could still be measured with the SIM telescope. The selec- 
tion 7 = was made because changing this value would 
result in a higher planetary mass and hence in a stronger 
astrometric signal, making the detection of orbital plane 
parameters even easier. Regardless of large error bars, we 
found that it is possible to detect the orbital plane parame- 
ters with Ta — 1-0 years. With this short timeline, the 99% 
error bars of the parameters were [-0.30, 0.22] and [0.39, 
1.54] for / and fi, respectively, demonstrating that it was 
indeed possible to determine their values. 

6. Conclusions and discussion 

The time needed to make a positive detection of an ex- 
trasolar planetary conpanion candidate depends essentially 
on its orbital period. It is commonly assumed that, to be 
able to detect the signature of such companion, an obser- 
vational timeline longer than the orbital period is required. 
Also, since most of the exoplanet candidates have been de- 
tected using the RV method, only the lower limit of their 
mass is available. With the aid of future space telescopes 
and accurate astrometric measurements, it will be possible 
to detect the inclination and thus the true mass of plane- 
tary candidates. 

We have shown that when high-precision RV and accu- 
rate astrometric measurements are both available, it is pos- 
sible to detect the true mass of stellar companions with ob- 
servational timelines considerably shorter than their orbital 
periods. Also, when the RV measurements have a long time 
span, astrometric measurements can reveal the true mass 
of a stellar companion in less time than one tenth of the or- 
bital period of the system. This ability is also demonstrated 
using the RV measurements of HD 154345 as an example. 
We find that, having these measurements with Try = 10.4 
years in hand, astrometric observations with SIM telescope 
are sufficient for obtaining the true mass, within a single 
year. 

Bayesian inference plays an important role when ex- 
tracting information from several sources of measurements. 



The ability to use RV and astrometric measurements simul- 
taneously makes it possible to employ observational time- 
lines below the orbital ones and still be able to make posi- 
tive exoplanet detections, thus helping to extract the max- 
imum amount of information from measurements and in- 
creasing the time efficiency of observations. 

In a forthcoming study, we plan to study the inclusion 
of additional transit-photometry measurements to further 
tighten the parameter probability densities in transiting 
scenarios. Also, the approach used here should be extended 
to systems with two or more planetary companions. 
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